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We consider the fractional oscillator being a generalization of the conventional linear oscillator in 
the framework of fractional calculus. It is interpreted as an ensemble average of ordinary harmonic 
oscillators governed by stochastic time arrow. The intrinsic absorption of the fractional oscillator 
results from the full contribution of the harmonic oscillators' ensemble: these oscillators differs a 
little from each other in frequency so that each response is compensated by an antiphase response 
of another harmonic oscillator. This allows to draw a parallel in the dispersion analysis for the 
media described by the fractional oscillator and the ensemble of ordinary harmonic oscillators with 
damping. The features of this analysis are discussed. 

PACS numbers: 05.40.-a, 05.60.-k, 05.40.Fb 



I. INTRODUCTION 

The harmonic oscillator, given by a linear differential 
equation of second order with constant coefficients, is a 
cornerstone of the classical mechanics (see, for example, 
[E0)- Today this elementary (and fundamental) concep- 
tion has the widest origin of physical, chemical, engineer- 
ing applications and needs no introduction. Its success 
mainly rests on its universality, and its simplicity gives 
boundless intrinsic capabilities for sweeping generaliza- 
tion. Suffice it to recall the passage from the language 
of functions in phase space to operators in Hilbert space 
so that the oscillatory model came strongly in the quan- 
tum theory [1, 0|- Therefore no wonder, the fractional 
calculus has made also an important contribution to this 
way. 

At first the approach had a formal character by chang- 
ing the second derivative in the harmonic oscillator equa- 
tion on the derivative of an arbitrary order. After find- 
ing out the solutions of such equations their relaxation- 
oscillation behavior was established [fl Q . The next step 
was a consideration of the total energy and the phase 
plane representation for the fractional oscillator [7|. To 
save the dimension of energy, it is necessary to generalize 
to the notation of momentum, though then the parameter 
m loses also the ordinary dimension of mass [8| . In this 
case the momentum is expressed in terms of the Caputo- 
type fractional derivative @. The fractional oscillator 
is like a harmonic oscillator subject to a damping. The 
source of the intrinsic damping is very intriguing. It is 
not evident from fractional calculus, from the generaliza- 
tion of derivative. The question requires an additional 
study exceeding the bounds of fractional calculus itself. 

Since it is a matter of the fractional integral/derivative 
with respect to time, the answer to the aforementioned 
problem should be sought by way of their concrete in- 
terpretation. Recently, the probability interpretation of 
the temporal fractional integral/derivative was suggested 



in [9(. There exists a direct connection between stable 
distributions in probability theory and the fractional cal- 
culus. The occurrence of the temporal fractional deriva- 
tive (or integral) in kinetic equations indicates the sub- 
ordinated stochastic processes. Their directional process 
is related to a stochastic process with a stable distribu- 
tion. The parameter characterizing the stable distribu- 
tion coincides with the index of the temporal fractional 
integral/derivative in the corresponding kinetic equation. 
This means that such a equation describes the evolution 
of a physical system whose time degree of freedom be- 
comes stochastic [l(|. The purpose of this paper is to 
expand the interpretation on the fractional oscillator. 

The paper is organized as follows. In Sec. |TI]we analyze 
an ensemble of harmonic oscillators with the stochastic 
time clock. The new clock (random process) substitutes 
for the deterministic time clock of the ordinary harmonic 
oscillator. The nondecrcasing random process arises from 
a self-similar a-stable random process of temporal steps. 
Using properties of the stochastic time clock, we obtain 
the equation for the fractional oscillator. In the spirit of 
this approach the fractional oscillator can be considered 
as an ensemble average of oscillators. Sec. Mil is devoted 
to the comparison of dispersion properties of the two me- 
dia. One of them consists of damped noninteracting har- 
monic oscillators, whereas another is the fractional oscil- 
lator. It turns out that their dispersion characteristics 
have a lot of common features. We discuss them in de- 
tail. Our conclusions are briefly summarized in Sec. IIVI 
Appendix contains calculations for the response of the 
driven fractional oscillator. They are useful for the dis- 
persion analysis in Sec. IIIII 



II. NORMAL MODES 

We start our consideration with the classical case of 
harmonic oscillator. Based on the Hamilton function 
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H = {p 2 
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''q 2 )/2, where p and q are the momentum 



and the coordinate, respectively, and to the proper fre- 
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quency, the motion equations take the form 

dq/dt = dn/dp = p, dp/dt = - dH/dq = -u 2 . (1) 

To multiply the first equation of ([1]) on ±iuj and to add 
it with the second equation, we arrive at 



The direct calculations give 



dc/dt = —iiuc, dc* / dt = iujc* 



(2) 



where the complex-conjugate values c and c* satisfy the 
relations 



c = (cuq + ip)/y2uj, c* = (ujq — ip)/y2to. 
The solutions of Eqs.® can be written as 



c(t) = c(0)e" 



[u>q(p)+ip(p)]e~**, 



c*(t) = c*(0)e-< = -=[uq(P) - ip(0)]e^. (3) 

VZUJ 



The values c and c* are called else the normal modes of 
oscillator . They have a very pictorial presentation in 
the form of a vector rotating just as, the hand revolves 
around the clock-face center with the frequency to. 

A physical system of harmonic oscillators coupled to an 
environment will interact with the environmental degrees 
of freedom. This leads to a damping of oscillatory mo- 
tion. If the interaction manifests itself at random fash- 
ion, one of possible ways to account for perturbations 
induced by the environment may be as following. Let 
us randomize the time clock of the value c(r) so that 
any characteristic time is absent. Assume that the time 
variable is a sum of random temporal intervals Tj on the 
non-negative semi-axis. If they are independent iden- 
tically distributed variables belonging to the strict do- 
main of attraction of a a-stable distribution (0 < a < 1), 
their sum has asymptotically (the number of the intervals 
tends to infinity) the stable distribution with the index 
a. Following the arguments of [T3, El, a new time clock 
is defined as a continuous limit of the discrete counting 
process N t = max{n g N : — where N is the 

set of natural numbers. The time clock becomes the hit- 
ting time process S(t). Its basic properties is represented 
in [ll|, [l2| ■ The probability density of the process S (t) is 
written in the form 



1 



Jilt — TU n ,Oi— 1 



u a - L du, 



(4) 



Br 



where Br denotes the Bromwich path. This probability 
density has a clear physical sense. It describes the prob- 
ability to be at the internal time r on the real time t. In 
the case we determine new normal modes 



c a (t) 



p S (t,r) c(r)dr, 



<&(*) = / P b (t,r) c*(r)dT. 
Jo 



Qa(t) 
Pa(t) 

where 
A(t) 

Bit) - 
and 



p(0) 



= [c* (t) + c a (t)]/V2ij = q(0)A(t) + ^B(t) , 

UJ 

= i[c* a {t) - c a {t)]^J/2= p(0)A(t) + wg(0)B(t) . 



p S {t,T) COS LOT d,T = E 2a ,l(-UJ 2 t 2a ) 



/"DO 

/ p s (t,r) smwrdr = ut a E 2a . a +i(~UJ 2 t 2a ) , 
Jo 



E, 



u^~ v du 
(u^ - z) 



is the two-parameter Mittag-Leffler function [i3j . Here 
it is easy to recognize the classical solutions for a = 1/2 
(exponential function) and a = 1 (sine and cosine). The 
functions A(t) and B(t) exhibit clearly the relaxation fea- 
tures for < a < 1/2, whereas for 1/2 < a < 1 the 
functions represent a damping oscillatory motion. The 
latter case just corresponds to the fractional oscillator. 
In particular the value A(t) satisfies the equation 



A(t) = A(0) 



T(2a) Jo 



(t-t') A{t')dt' , 



with A(0) = 1, where T(x) denotes the gamma function. 
The appropriate equation can be written also for B(t). 
It should be recalled here that the power kernel of frac- 
tional integral of the order a, < a < 1 "interpolates" 
the memory function between the Dirac <5-function (the 
absence of memory) and step function (complete ideal 
memory). This means that such memory manifests itself 
within all the time interval (0, t), but not at each point of 
time (complete but not ideal memory). Under the ideal 
complete memory the system "remembers" all its states, 
and this excites the harmonic oscillations in such a sys- 
tem. The absence of memory causes only the relaxation. 
The order of fractional integral represents a quantitative 
measure of memory effects [14j . In accordance with the 
theory of memory effects the fractional oscillator contains 
simultaneously the oscillatory motion and the relaxation. 

From the series representation of E^ v (z) we derive 
the leading asymptotic behavior of the value A(t) and 
B(t) for t -> 0: \\m t ^ A{t) = 1, lim t _y B{t) = 0. Ac- 
cording to [13], the two-parameter Mittag-Lcfflcr func- 
tion approaches zero as z — > oo in the sector of angles 
| arg(— z)\ < (1 — fi/2)n, and increases indefinitely as 
z — > oo outside of this sector. In our case we can use 
the following expansion valid on the real negative axis 



JV-l 



E lhV (z) = -Y, 



r(i/ - nfi) 



+ 0(1*1""), z 



Thus, for < a < 1/2 and 1/2 < a < 1 the value 
A{t) and B(t) decreases algebraically in time. As dis- 
tinct from the case of a damping harmonic oscillator, the 
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model describes another damping mechanism, without 
any external frictional force. The damping of a fractional 
oscillator is due to internal causes [15| . How to explain 
the attenuated oscillations? This important feature of 
fractional oscillator has been already noted from time to 
time in various publications 0,0)9- However, the source 
of such intrinsic damping remained undecided. 

We suggest the following interpretation. The fractional 
oscillator should be considered as an ensemble average 
of harmonic oscillators. When all harmonic oscillators 
are identical, and we set their going in the same phase, 
their full contribution will be equal to the product of 
the number of oscillators and the response of one oscil- 
lator. This occasion appears if a = 1. However, if the 
oscillators differ a little from each other in frequency, 
even if they start in phase, after a while the oscillators 
are allocated uniformly up to the clock-face. Each re- 
sponse will have an antiphase response of another oscil- 
lator so that the total response of all harmonic oscillators 
in such a system is compensated. Although each oscil- 
lator is conservative (its total energy saves), the system 
of such oscillators, resulting in the fractional oscillator 
(1/2 < a < 1), shows a dissipative nature. In this con- 
nection it should be pointed out that the similar situation 
may be observed also in the medium of harmonic oscilla- 
tors, having a given probability density on frequency (for 
example, the Lorentz distribution [Hj]). Both these cases 
are closely connected with each other and have a com- 
mon ground, though, generally speaking, they describe 
different physical systems. As has been shown in [l7lll8j. 
Lagrangian and Hamiltonian mechanics formulated with 
fractional derivatives in time can be used for the descrip- 
tion of nonconservative forces such as friction. It should 
be also mentioned an interpretation of fractional oscilla- 
tor in [19]. In this case the Liouville equation is formu- 
lated from a fractional analog of the normalization con- 
dition for distribution function that can be considered 
in a fractional phase space. The latter has a fractional 
dimension as well as the fractional measure. The vol- 
ume element of the fractional phase space is realized by 
fractional exterior derivatives. The usual nondissipative 
systems become dissipative in the fractional phase space. 
However, the approach is different on ours. It operates 
with fractional powers of coordinates and momenta. Such 
fractional systems are nonlinear. 



and x(0) = is described by the following equation 



x ® = "rTT / (t-tT- 1 x(t')dt> + 
r(a) J 



i r* 



r(a) Jo 



{t-ty- 1 F{t')dt' 



(5) 



where it should be kept 1 < a < 2, and F is the exter- 
nal force. The dynamic response of the driven fractional 
oscillator was investigated in [§[: 



x(t) 



F(t') (t - t')*- 1 E a , a (-u>g(t - t') a ) dt' . (6) 



This allows us to define the response for any desired forc- 
ing function F(t). The "free" and "forced" oscillations 
of a fractional oscillator depend on the index a. How- 
ever, in the first case the damping is characterized only 
by the "natural frequency" ujq, whereas the damping in 
the case of "forced" oscillations depends also on the driv- 
ing frequency u>. Each of these cases has a characteristic 
algebraic tail itself, associated with damping 15]. 

Let F(t) be periodic, F(t) = F a e Jwt . Then the solution 
of Eq. ([5]) is determined by taking the inverse Laplace 
transform 



x(t) 



1 

2^7 



Br 



Fp (s + juj) ds 
(s 2 + io 2 ) (s a + wg<) 



(J) 



The Bromwich integral ([7]) can be evaluated in terms of 
the theory of complex variables. Some particular exam- 
ples of forcing functions were considered in Q. However, 
the set turns out to be enough scanty for our aim. The 
necessary computations with F(t) = Asm(ujt + <j>) are 
fulfilled in Appendix. The phase <f> is constant. 

If one waits for a long enough time, the normal mode 
of this system is damped. Therefore, consider only the 
forced oscillation. After the substitution of x(t) — Xo< 
for x(t) in Eq.© we obtain 



Jut 



x e 



jut 



III. DISPERSION 



"(a) Jo 

+ tJ~t [ (t-tT-'Foe^'dt'. (8) 
r ( a ) Jo 



It is convenient to change the variable u>(t — t') = £ in the 
integrand. Next we can divide out exp(jwt) from each 
side of (|8]) and direct t to infinity. The procedure permits 
to extract the contribution of steady-state oscillations. 
Using the table integral 2Q\ 



Now we examine the behavior of the fractional os- 
cillator under the influence of an external force. From 
above this case corresponds to oscillations in the ensem- 
ble of nonidentical harmonic oscillators noninteracting 
with each other. In the framework of this model the 
fractional oscillator with the initial conditions x(0) = 



z a - x e^ z dz = T(a)e-^ a/2 



Eq. 



gives 



Fn 



Xo 



K + u a exp(j7ra/2)] 



(9) 
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FIG. 1: Dispersion dependence of the fractional oscillator in 
the form of the functions fF(uj/oJo,o) and qf{oj/ujo, a) with 
different values a, from 0.1 to 0.9 with the step 0.1. 



FIG. 2: Dispersion dependence in the classical case (ensemble 
of ordinary harmonic oscillators with damping) with different 
values /3, from 0.1 to 1.0 with the step 0.1. 



This result is completely confirmed by a more compre- 
hensive analysis given in Appendix. 

As is well know, the ensemble behavior of identical 
noninteracting harmonic oscillators is a basic topic for 
considering in the classical theory of dispersion. The non- 
identity of oscillators is necessary to take into account, for 
example, for the dispersion analysis of propagating elec- 
tromagnetic waves into a heated gas, where the spread 
in molecule velocity values leads to a Dopplcr shift of the 
oscillators' normal frequency with respect to the forced 
held frequency. Right now let a medium of oscillators be 
such that results in the fractional oscillator. It interests 
us the polarizability of such a medium. In this case the 
permittivity is written as e = 1 + 4.7re 2 :co/(-Fo TO ), where 
e is the electron charge. It should be pointed out that 
in contrast to a simple harmonic oscillator the param- 
eter to does not have the ordinary dimension of mass. 
However, the generalized momentum p of the fractional 
oscillator is defined via the Caputo-type fractional deriva- 
tive of order a/2 @ so that the expression p 2 /(2m) has 
the dimension of energy (see details in Q). The real and 
imaginary parts of the permittivity take the form 



Ree(w) = 1 H 
Ime(w) = 



4ne 2 [uj% + u a cos(7ra/2)] 
to [io 2a + Lu 2a + 2u%uj a cos(?ra/2); 

47re 2 w a sin(7ra/2) 
" to [uj 2a + uj 2a + 2wJw Q cos(vra/2)] 



,(10) 



(11) 



For a = 2 we arrive at the Sellmeier's formula 12111 



e(uj) = 1 + 4ire 2 N/[m{u) 2 



u 2 )] 



where we include N to account for the number of har- 
monic oscillators in the medium. In this case the pa- 
rameter to is really the electron mass. The index a = 2 
corresponds to the classical harmonic oscillator without 



any damping, and all the oscillators in the ensemble go 
in the same phase. Therefore, the Sellmeier's formula 
contains only the real part of permittivity. 

We can conduct a clear comparison between the disper- 
sion characteristics of the fractional oscillator and one of 
an ensemble of classical harmonic oscillators with damp- 
ing. Normalize the frequency lo in their permittivity 
by loq. In fact the constants (like e, to and so on) in 
Ree(w) and Ime(u;) define only a scale. Thus, one can 
pick out the functional dependence of these permittivi- 
ties on lu/ujq = z. Denote 2j/luq by /?, where 7 defines 
the damping in each classical harmonic oscillator. Then 
we have the following dependences for the fractional os- 
cillator 



1 + z a cos(7ra/2) 
z 2a + 2z a cos(7ra/2) + 1 ' 

z a sin(7ra/2) 
z 2a + 2z a cos(vra/2) + 1 ' 

and ones for the classical harmonic oscillators with damp- 
ing 



HeeF(uj) — > fp(z,a) 
lme F (oj) -t gF{z,a) 



Ree D (uj) 
Im c d (uj) 



f D {z,(3) 



l-z 2 



(1- z 2 ) 2 + z 2 f3 2 
zp 

(l-z 2 ) 2 + z 2 /3 2 



If the parameter j3 determines the damping value in the 
harmonic oscillator, the index a just characterizes the 
same for the fractional oscillator. The extremum values 
of Jd{z, p) and <7n(z, /?) decrease with increasing the pa- 
rameter /3 whereas for fp(z,a) and gp(z,a) vice versa: 
the extremum values increase with increasing the index 
a, though it should be noted that this index itself belongs 
only to the interval 1 < a < 2. The functions f...(z) and 
g...(z) are shown on Fig. [5] and Fig. [3J 
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From the relations (fTt)|) and (fTTj) it follows that there is 
a frequency range, where the absorption is small, and the 
refraction coefficient increases with frequency (normal 
dispersion). Moreover, in the frequency range, where the 
absorption is big, the anomalous dispersion happens to 
be the case for the refraction coefficient decreasing with 
frequency. In this connection it should be pointed out 
that the presence of the normal and anomalous disper- 
sion is typical for such an ensemble of ordinary harmonic 
oscillators and is well known. However a new fact estab- 
lished here is that the normal and anomalous dispersion 
is also typical for the medium described as a fractional 
oscillator. 



IV. SUMMARY 

We have shown that the fractional oscillator can 
be considered as a model of the harmonic oscillators' 
medium. Its stochastic properties accumulate in the in- 
dex of the fractional integral/derivative with respect to 
time. The frequency distinction of the oscillators (con- 
stituents of the fractional oscillator) from each other is 
at the bottom of the intrinsic damping for such a sys- 
tem. As a consequence, the dispersion properties of the 
medium, like the fractional oscillator, is enough similar 
to the case, when a medium is modeled by an ensemble 
of harmonic oscillators with damping. 




FIG. 3: Contour, inside which the function x(s) remains 
single-valued and analytical all over, with the exception of 
poles and ljq exp(±jir /a). 



Appendix 

We here derive properties of the response function (|6]) 
for the forcing function Asm(ut + 4>) directly from its 
representation as a Laplace inverse integral 

x(t) = I e st x(s)ds = 

2tt.7 J Br 

1 f st A (s ■ sin0 + u> ■ cos</>) ds 



2nj 



Br 



( S 2 + OJ 2 ) {S a + W%) 



, (Al) 



where the phase (j) is constant, Br denotes the Bromwich 
path, and 1 < a < 2. By bending the Bromwich path 
into the equivalent Hankel path (Fig. [3]), the response 
function x{t) can be decomposed into two contributions. 

The first contribution arises from the two borders of 
the cut negative real axis (lines DE and FG): 

if i r°° 

xi(t) = -- — / e st x(s)ds-- — / e st x(s)ds. 
27T? J_ 0O 2nj J 

To enter s = re J7r into the integral taken along the upper 
border and s — re~ J ™ into the integral along the lower 
border, we get 



x x (t) = / e- rt M a {r)dr 
Jo 



with 

M a (r) 



A r a (w cos 4> — r sin (ft) sin ira 
it (r 2 + U! 2 ) (r 2a + 2r a wg cos na + w§ Q ) ' 



The second contribution is determined by the Cauchy 
theorem on residues. The integrand of (Al) has the fol- 
lowing poles 

s = ±ju> and s = woe ± -' n ^ a . 

Calculating the residues of the poles s = ±ja;, we obtain 

, A sin(wt + 0) + w Q sin(wi — 7ra/2 + (/)) 
X2 ^' ~ [ u 2a + u) 2a + 2uj*uj a cos(7ra/2) _' 

It remains to define the residues for the two other poles: 
e st (s ■ sin0 + oj ■ cosc/>) 



a ujQt{cos 7r /aij sin 7r/a) 



w e 



±jT/a s i n ^)-f wcos< 



IW " 1 e ±jV(a-l)/a [cj2 e ±2j7r/a + 



They lead to 



4(t) 



2A exp(woi cos(7r/a)) [Ccos0 — D sin <p] 
aw a - 1 (^+w 4 + 2a;^ 2 cos(27r/Q!)) ' 



6 



where 

C = cj[cjq cos(cj t sin(7r/a) — 7r(l + a)/a) + 

+ lu 2 cos(u;oi sin(-7r/a) + n(l — a)/a)], 

D = o;o[o;q cos(a;ot sin(7r/a)) + 

+ lu 2 cos(woi sin(7r /a) + 2tt /a)] . 

As a result, the response function x(t) takes the form: 

x(t) = xi(t) + x' 2 (t) + x' 2 \t) . 

Since cos(7r/a) < 0, the term x 2 (t) describes the relax- 
ation of the normal mode in this system. For 1 < a < 2 
and all r the denominator of the value M a (r, a) is always 
positive: (r 2a + 2r a uig costfc* + ujI") > (r a - ujg) 2 > 0, 
and the term sin7ra is always negative. Depending on cj>, 
each of terms ui cos <f> and r sin (/> may be both positive 



and negative. However the value x\(t) becomes vanish- 
ingly small with t — > oo. The steady-state oscillation in 
this system is defined only by the term x' 2 {t). The latter 
can be expressed as x' 2 (t) = A\ sin(wt + — 5), where 



[uj 2a + uj 2a + 2uj a u% cos(7ra/2)] 1 /2 ' 

[ u a sin(7ra/2) 

o = arctan ; — . 

_L0 a cos(7ra/2) + uj% _ 

To put 4> = in (Al), we arrive at the results of section 
4.3 from ||. It should be also noted that the oscilla- 
tory contribution X2(t)\<f,=o has some resemblance with 
the "free" oscillations of a damped harmonic oscillator 
and the forced oscillations of a driven damped harmonic 
oscillator [l5j . 
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